Spin dynamics in molecular ring nanomagnets: Significant effect of acoustic phonons 

and magnetic anisotropics 
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The nuclear spin-lattice relaxation rate 1 /Ti is calculated for magnetic ring clusters by fully diag- 
onalizing their microscopic spin Hamiltonians. Whether the nearest-neighbor exchange interaction 
J is ferromagnetic or antiferromagnetic, 1/Ti versus temperature T in ring nanomagnets may be 
peaked at ksT ~ |J| provided the lifetime broadening of discrete energy levels is in proportion to 
T 3 . Experimental findings for ferromagnetic and antiferromagnetic Cu 11 rings are reproduced with 
crucial contributions of magnetic anisotropies as well as acoustic phonons. 

PACS numbers: 75.50.Xx, 75.40.Mg, 75.75.+a, 76.60. -k 



In recent years, considerable efforts pj have been de- 
voted to constructing and investigating magnetic systems 
of nanoscale dimension that comprise a controllable num- 
ber of transition metal ions. Resonant magnetization 
tunneling 0, in [Mn^O^CHsCOO^HaO^], now 
well-known as Mnl2, has stimulated increasing interest 
in mesoscopic magnetism. Among such topical molecu- 
lar nanomagnets are highly symmetrical clusters of al- 
most planar ring shape 0, 0| . They are highly varied, 
containing four to eighteen metal ions of spin ^ to | and 
therefore, serve to reveal a quantum-to-classical crossover 
between molecular and bulk magnets [||. 

Nuclear magnetic resonance (NMR) is an effective 
probe of low-energy spin dynamics, and the nuclear spin- 
lattice relaxation time Ti ha s been measured for various 
molecular wheels II I1M 13 13 13 13 El Interest- 
ingly, 1/Ti as a function of temperature T under a fixed 
field H is commonly peaked at k-gT ~ |J|, where J is 
the intracluster exchange interaction between neighbor- 
ing ion spins. Baek et al. ^3 have recently provided a 
key to this long-standing problem by reanalyzing exten- 
sive observations of antiferromagnetic rings. When 1/Ti 
is divided by the static susceptibility-temperature prod- 
uct xT, its peak is more pronounced and well-fitted to 
the Lorentzian-type expression, uj c (T)/[uj^(T) +u^(H)}, 
where = "f^H is the Larmor frequency of probe nu- 
clei, whereas ui c is what they define as the temperature- 
dependent correlation frequency. The renormalized re- 
laxation rate is thus peaked at a temperature satisfying 
uj c (T) ~ uj-n(H). Considering the significant difference 
between the electronic and nuclear energy scales (Tiwn 
10~ 5 |J|), huj c (T) may be ascribed to the averaged life- 
time broadening of discrete energy levels. Demonstrat- 
ing that chromic and ferric wheels give uj c oc T a with 
a = 3.0 ~ 3.5, Baek et al. claim that the exchange- 
coupled ion spins are likely to interact with the host 
molecular crystal through the Debye-type phonons. 

In response to this stimulative report, several authors 
[lil H?| inquired further into the underlying scenario. 
While their arguments were elaborately based on micro- 
scopic spin Hamiltonians and enlighteningly verified the 



relevance of spin-phonon coupling to the notably peaked 
1/Ti, any magnetic anisotropy was neglected and/or 
most of the transition matrix elements were discarded 
in their evaluation of the dynamic spin correlation func- 
tions. Thus, we take the naivest but thus cumbersome 
approach: We set up realistically dressed Hamiltonians, 
completely diagonalize them, and sum up all the tran- 
sition matrix elements into the dynamic structure fac- 
tor. Such calculations are inevitably restricted to suffi- 
ciently small clusters but can nevertheless illuminate key 
factors in nanoscale spin dynamics, intrinsic intracluster 
anisotropies and extrinsic intercluster phonons. 

We consider both ferromagnetic and antiferromagnetic 
rings of various metal ion spins and describe them by the 
Hamiltonian 



H = J^2s r S l+1 +J2d r (S l x S l+1 ) 
i=i i=i 
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where Sl+i = The first term is the isotropic Heisen- 
berg exchange interaction. The second term introduces 
the antisymmetric Dzyaloshinsky-Moriya (DM) interac- 
tion with di set for (— l) l (d, 0, 0). The third and fourth 
terms describe the axial magnetic anisotropy in the cases 
of ferromagnetic (J < 0) and antiferromagnetic (J > 0) 
exchange coupling, respectively. Both Z?fm and -Dafm 
are simply written as D in all the figures. The last term 
represents the Zeeman interaction with an external field. 
The nuclear spin-lattice relaxation rate is generally ex- 
pressed as 
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\Aq T \ 2 is the form factor describing the hyperfine cou- 
pling of spin excitations at wave vector q with the probe 
nuclei. S aT (q, to) is the dynamic spin structure factor and 
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is given by 

S"{q,u;) =Y^e- n ^ k - T (i\S^\j){j\S%\i) 

x&iu + Ui -w i )/X) e ~' iW ' /fcBT ' ( 3 ) 

where S° is the Fourier transform of Sf and \i) is an exact 
eigenstate of the Hamiltonian (JTJ with energy fto;*. There 
may be longitudinal contributions S zz (q, ojn) to 1/Ti in 
general, originating from anisotropic hypcrfinc coupling, 
and therefore strongly depending on the detailed molec- 
ular structure. Considering that the characteristic peak 
of 1/Ti versus T is commonly observed for various ring 
nanomagnets, here, we assume that the transverse spin 
fluctuations S~* (q,u>Tss) predominate in the relaxation 
rate ®. We further replace A+~ by A+T Q = A + ~ [ftllll|. 
because the geometric coefficients A^ T usually exhibit lit- 
tle momentum dep endence, particularly in proton NMR 
measurements |l6l Il8j . In the limit of H — ► Q without any 
anisotropy, (q,w) is equivalent to S zz {q,ui) and can 
be obtained rather easily. However, it is not the case with 
any NMR measurement. Even in the case of d = 0, where 
the Hamiltonian Q commutes with the total magnetiza- 
tion and a translation by one site along the ring, it takes 
a long time to straightforwardly evaluate (q, u>) for a 
hexanuclear Fe m ring (Fe6) or an octanuclear Cr m ring 
(Cr8). With nonvanishing DM interaction, the naivest 
treatment of Fe6 and Cr8 with the Hamiltonian is 
impossible to begin with. Therefore, we try to extract 
the essential physics from more tractable clusters. 

The most important ingredient in calculating eq. (0) 
is a treatment for 5(u> + LOi — u>j). When we adopt the 
Lorentzian-type expression 
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hr may be recognized as the lifetime broadening of the 
discrete energy levels ^^|- There has been a pioneering 
attempt |l7| to solve a coupled spin-phonon Hamiltonian, 
leading to a level-dependent expression of T. However, 
such a calculation artificially selects particular transition 
matrix elements out of the summation in eq. © and is 
feasible in a rather limited situation, for instance, at suf- 
ficiently low temperatures &bT <§C Swd, where Swo/fe is 
the Debye temperature, which is usually less than 100 K 
for molecular nanomagnets 0|. Here, we place great 
emphasis on including full spin degrees of freedom into 
our calculation, and therefore discuss T in an empiri- 
cal way. We set T proportional to T 3 on the one hand 
and independent of temperature on the other hand. The 
former treatment is motivated on the basis of the exper- 
imental observations by Baek et al. 0] and results in 
supporting their scenario — discrete energy levels of anti- 
ferromagnetic rings are lifetime-broadened owing to spin- 
acoustic phonon coupling. The latter treatment attempts 
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FIG. 1: (Color online) 1/Ti versus T with varying fields ap- 
plied to Cu6 (a), Mn4 (b), Cu8 (c), and Ni8 (d), which are 
assumed to be free from any anisotropy. Red (toned-down) 
and black lines show calculations on the assumption that T 
is independent of T and proportional to T 3 , respectively. 
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FIG. 2: (a) Temperature at which 1/Ti reaches its maximum 
is plotted versus H. (b) Maximum value of 1/Ti is plotted 
versus 1/H, where any anisotropy is set equal to zero, and T 
is assumed to be proportional to T 3 . 



to exclude any possibility of phonons assisting the nuclear 
spin-lattice relaxation and results in suggesting more in- 
trinsic mechanism such as intercluster dipolar interac- 
tions [2(| for the lifetime broadening of discrete energy 
levels. 

First, we consider various magnetic rings without any 
anisotropy in an attempt to elucidate the effects of acous- 
tic phonons on the relaxation rate. Figure ^ shows 
the temperature dependences of 1/Ti f° r ferromagnetic 
hexanuclear Cu 11 (Cu6) plf. ferromagnetic tetranuclear 
Mn m (Mn4) [2^, antiferromagnetic octanuclear Cu 11 
(Cu 8) |2jj, and antiferromagnetic octanuclear Ni 11 (Ni8) 
[24| rings. When we set T to be independent of T, 
as well as H, the resultant 1/Ti is reminiscent of \T; 
that is to say, it monotonically decreases and increases 
with T in the ferromagnetic and antiferromagnetic cases, 
respectively. No peak appears in the curve of 1/Ti 
versus T except for a trivial ceiling due to the Zee- 
man splitting of the ferromagnetic ground state. If 
we adopt T to be proportional to T 3 , a notable peak 
appears at T = (T) pca k ~ |J|/fce, which is qualita- 



3 



Cu6: Tccr-' 



(a)^ B ///l7l = 3.64xlO- 3 

D/J = 

D/J =0.003 

\ D/J= 0.007 

— D/y =0.012 




(b)^ B ///UI = 7.35xl0- J 

— o/y=o 

Z5/y =0.003 

D/J =0.007 

zj/y =0.012 



..(c) gfl B tf/UI= 1.62x10"- 




(d)g|l B tf/|y| = 3. 13x10^ 

<----. A — o/y=o 



■D/J = 0.003 
■B/7 = 0.007 
= 0.012 



g// B «/|y| = 3.64xl0- 3 
g/J B ///|y| = 7.35xlO- 3 
g// B ///|y|=1.62xlO- 2 
g/j B ///|y|=3.13xl0- 2 



(e)D/J= 0.007 



1 2 

fc B T/iyi 



1 2 
fe,7V|yl 



1 t B r/iyi 2 



fc R 77|yl 



1 2 

* B r/ui 



FIG. 3: 1/Ti versus T for Cu6 with anisotropic exchange interaction, where F is set proportional to T 3 
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FIG. 4: 1/Ti versus T for Cu8 with antisymmetric interaction, where F is set proportional to T 3 
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tively consistent with numerous experimental observa- 
tions 0, H 0, 113, EE EE 13 We have confirmed that 
such a peak structure does not appear with any other 
temperature dependence of r, whether it be of power- 
law type or of activated type. With increasing field, the 
peak becomes smaller and shifts to higher temperatures 
monotonically. Figure Ufa) shows (T) pea k versus H and 
reveals its power-law behavior (T) pca k oc H v . We find 
that v ~ 0.32, 0.34, 0.30 and 0.31 for Cu6, Mn4, Cu8 
and Ni8, respectively. Considering the significant differ- 
ence between the nuclear and electronic energy scales, 
Tk^N ~ 10 -5 J, it must be the averaged lifetime broaden- 
ing of discretized energy levels, referred to as TiWc [l2|. 
that characterizes 1/Ti. The relaxation should be most 
accelerated at = H ~ uj c ; that is, we suppose that 
1/Ti is proportional to lo c /{uj^ +uj^). Then, the above 
findings indicate that tu c oc T a , with a ranging from 2.9 
to 3.3, which is in very good agreement with some exper- 
imental findings |l2j . However, there are observations of 
(T)peak remaining almost unchanged with H H it 
is also interesting to plot the maximum value of 1/Ti, 
(1/Ti) max , versus 1/H. Figure [^b) shows that in the 
ferromagnetic cases, there well holds a linear relation be- 
tween them, which is again consistent with the above con- 
sideration, 1/Ti oc uj c /(ui^ -fw^) < 1/2-f^H, whereas in 
the antiferromagnetic cases, the calculations deviate from 
such a straightforward relationship. The peak structure 
indeed depends not only on the constituent metal ion 
spins, but also on their surrounding ligands 0. We are 
thus led to consider the effects of magnetic anisotropy. T 
is always set proportional to T 3 in the following. 

The ferromagnetic ring Cu6 is simulated in more de- 
tail in Fig. 21 where the anisotropic exchange interaction 
T>fm is taken into account. The anisotropy effect is so sig- 



nificant at fceT ~ | J [Figs. EI a ) — EJd)] as to completely 
break the monotonic field dependences of both (T) pea k 
and (1/Ti) max [Fig. EHe)]. Even a double peak may ap- 
pear in 1/Ti versus T. Metallic wheels jH El El Q 
can act as hosts for an alkali-metal ion, which affects the 
coordination geometry of the environmental transition- 
metal ions. The experimental findings for hexagonal fer- 
ric wheels with a lithium ion, a sodium ion and no alkali 
ion at the center are indeed different from each other in 
a subtle way El an d such observations are consistent 
with the sensitivity of 1 jT\ to magnetic anisotropy. 

The antiferromagnetic ring Cu8 is free from the crys- 
talline anisotropy Z?afm- On the other hand, 63 Cu nu- 
clear quadrupole resonance spectra of Cu8 rings indicate 
four crystallographically inequivalent Cu 11 ions in each 
ring, which is attributable to the slight deviation of Cu 11 
ions from a planar octagon J2j| and/or possible distortion 
of the octagonal environment due to spin-lattice coupling 
|29j. Hence, we introduce a DM term into the Hamil- 
tonian of Cu8 and show the resultant l/Ti in Fig. 0] 
The DM interaction induces another peak at a tempera- 
ture much lower than J/k-Q [Figs. EIa)-0Id)]. As soon 
as d is switched on, there appears a renewed peak at 
fcgT ~ 0.4J. With increasing d, the peak remains sit- 
ting at almost the same temperature, whereas its height 
monotonically increases and reaches saturation. Thus, a 
DM term in the Hamiltonian masks the peak structure at 
fceT ~ J. Both (T) pea k and (1/Ti) max for the renewed 
peak monotonically decrease with increasing field, but 
their field dependences are much more moderate than 
those of the original peak [Fig. 0I e )] ■ 

Finally, we attempt to reproduce experimental find- 
ings for the ferromagnetic hexanuclear Cu 11 rings 
[(PhSi0 2 ) 6 Cu 6 (0 2 SiPh) 6 ]-6EtOH H2 and the antiferro- 
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magnetic octanuclear Cu 11 rings [Cu8(C5H7N 2 )8(OH) 8 ] 
[23j within the framework of our theory. We consider an 
anisotropic exchange contribution D/\J\ = 0.007 with 
J/k-Q = — 60. 5 K for the Cu6 clusters [21|, while alter- 
nating DM vectors d/J = 0.015 with J/fc B = 800 K 
for the Cu8 clusters 0. Figures |3{ a) and EJb) show 
the thus-calculated 1/Ti for the Cu6 and Cu8 clus- 
ters, respectively, where the lifetime broadening of dis- 
crete energy levels is taken as F — c(/cbT /|J|) 3 with 
c = 3.2 x 10 10 rad • Hz (Cu6) and c = 1.0 x 10 11 rad • Hz 
(Cu8), while the hyperfme coupling strength is set for 
(7nA+-) 2 = 18 x 10 14 rad 2 Hz 2 (Cu6) and {j N A+-) 2 = 
2.8 x 10 14 rad 2 Hz 2 (Cu8). The present coupling constants 
are so reasonable as to fit with the experimental esti- 
mates .9]: (7nA+-) 2 = (8 ±2) x 10 14 rad 2 Hz 2 (*H NMR 
on Cu6) and ( 7N A+-) 2 = (9 ± 3) x 10 14 rad 2 Hz 2 (*H 
NMR on Cu8) within a factor of about three. In the Cu6 
clusters, 1/Ti is sensitive to the anisotropic exchange in- 
teraction at k^T < 2 1 J| in particular, and that is why 
we find such a reduced field dependence of (1/Ti) max 
in Fig. Eta). The ratio of (1/Ti) max (# = 0.164 T) to 
(1/Ti) max (i/ = 1.409 T) is two or less, whereas in an- 
other measurement on the same material @ the ratio 
is reported to reach three. Considering that both mea- 
surements were performed on polycrystalline powders, 
such a sample dependence of T\ findings sounds con- 
vincing. There exists a single crystal of Mn4 ferromag- 
netic rings instead [22^ . where any anisotropy effect can 
be verified quantitatively. In the Cu8 clusters, the DM 
interaction appears to be essential for peaks to appear 
at much lower temperatures than J/k^. The measure- 
ment is unfortunately limited to &bT < 0.5J, because 
the Cu8 sample breaks down at above room temperature. 
The observations nevertheless suggest that 1/Ti almost 
reaches saturation at room temperature. The probable 
low-temperature peaks in the Cu8 clusters may be char- 
acteristic of an antisymmetric interaction and should be 
distinguished from peaks originating from the exchange 
interactions, whether isotropic or anisotropic, which usu- 
ally appear at k^T ~ J. 

We have microscopically interpreted the temperature 
dependences of 1/Ti in molecular ring nanomagnets. In 
comparison with antiferromagnetic rings 0> EH EH EH 
E3, ferromagnetic ones have so far much less been cal- 
culated microscopically. Many of antiferromagnetic fea- 
tures were derived from isotropic Hamiltonians and were 
hardly compared with corresponding experimental ob- 
servations. In such circumstances, we have investigated 
both ferromagnetic and antiferromagnetic rings with par- 
ticular emphasis on possible anisotropic terms in their 
Hamiltonians and their full diagonalization. While 1/T\ 
versus T in exchange-coupled ion spins of ring shape 
may be notably peaked at k^T ~ |J| provided they 
are coupled to acoustic phonons of the host molecular 
crystal, the peak structure is very sensitive to magnetic 
anisotropics. With increasing anisotropies due to single- 
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FIG. 5: (Color online) 1/Ti versus T for proton 
nuclei in [(PhSi0 2 )6Cu 6 (02SiPh) 6 ] ■ 6EtOH (aj_ and 
[Cu 8 (C5H7N 2 )8(OH) 8 ] (b). Experimental findings U EH 
(symbols) are compared with calculations with (black lines) 
and without [red (toned-down) lines] anisotropies, where F is 
set proportional to T 3 . 



ion and dipolar interactions, which are all describable 
within symmetric spin-bilinear terms, both (1/Ti) max 
and (T) pca k move in an irregular way. With increas- 
ing antisymmetric DM interaction, on the other hand, 
another peak grows monotonically far from the original 
peak at k^T ~ J. The sample- dependent field depen- 
dence of 1/Ti for Cu6 clusters |3, III and the almost 
saturated 1/Ti at room temperature for Cu8 clusters 
0, ^3 are both understandable in this context. Indeed 
anisotropic contributions to the Hamiltonian are no more 
than a few percent of the exchange coupling constant in 
most cases |2lL I33I I34I EH , but they potentially cause a 
drastic change in 1/Ti particularly at low temperatures 
fceT ^ I J\. With increasing field, the 1/Ti peak becomes 
smaller and shifts to higher temperatures in Cr8 antifer- 
romagnetic rings |l2j, which is reminiscent of Figs. [T^c) 
and IHd) , while the peak remains sitting at almost the 
same temperature in Fe6 antiferromagnetic rings El Bl , 
which may be ascribed to alternating DM vectors |l3u34| . 
With nonvanishing DM vectors, neither Cr8 nor Fe6 clus- 
ters can be fully diagonalized. Low-temperature observa- 
tions, such as quantum fluctuations of the ground-state 
spin multiplet |3j| and level crossings in an applied field 
[37J, may be understandable within an effective theory 
based on selected low-lying states, whereas the essentially 
thermal behavior of our interest should be examined with 
the whole Hamiltonian. Understanding will increase with 
further measurements of [MruCUfCgHgNO^ ^I |22^ and 
[(C 2 4H 16 N 8 2 )4Ni 8 (H 2 0) 8 ](BF4)8 • 16H 2 j3, which 
are ideal ferromagnetic and antiferromagnetic rings, re- 
spectively, with relatively small spin degrees of freedom 
and moderate intramolecular exchange interactions. 

We thank S. Miyashita, Y. Furukawa, F. Borsa, and S. 
Maegawa for fruitful discussions. We give our warmest 
acknowledgments to Y. Furukawa and F. Borsa for pro- 
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ported by the Ministry of Education, Culture, Sports, 
Science, and Technology of Japan. 
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